CERN-TH/2003-073 
|hep-ph/0303242| 



Sneutrino Inflation in the Light of WMAP: 
Reheating, Leptogenesis and Flavour- Violating Lepton Decays 



John Ellis 1 , Martti Raidal 12 and T. Yanagida 3 

1 TH Division, CERN, CH-1211 Geneva 23, Switzerland 
2 National Institute of Chemical Physics and Biophysics, Tallinn 10143, Estonia 
3 Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 



ABSTRACT 



We reconsider the possibility that inflation was driven by a sneutrino - the scalar supersym- 
metric partner of a heavy singlet neutrino - in the minimal seesaw model of neutrino masses. 
We show that this model is consistent with data on the cosmic microwave background (CMB), 
including those from the WMAP satellite. We derive and implement the CMB constraints on 
sneutrino properties, calculate reheating and the cosmological baryon asymmetry arising via 
direct leptogenesis from sneutrino decays following sneutrino inflation, and relate them to light 
neutrino masses. We show that this scenario is compatible with a low reheating temperature 
that avoids the gravitino problem, and calculate its predictions for flavour-violating decays of 
charged leptons. We find that \i — > should occur close to the present experimental upper 
limits, as might also r — > /ry. 
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1 Introduction 



Inflation pQ has become the paradigm for early cosmology, particularly following the recent 
spectacular CMB data from the WMAP satellite 2], which strengthen the case made for infla- 
tion by earlier data, by measuring an almost scale-free spectrum of Gaussian adiabatic density 
fluctuations exhibiting power and polarization on super-horizon scales, just as predicted by 
simple field-theoretical models of inflation. As we review below, the scale of the vacuum energy 
during inflation was apparently ~ 10 16 GeV, comparable to the expected GUT scale, so CMB 
measurements offer us a direct window on ultra-high-energy physics. 

Ever since inflation was proposed, it has been a puzzle how to integrate it with ideas in particle 
physics. For example, a naive GUT Higgs field would give excessive density perturbations, and 
no convincing concrete string-theoretical model has yet emerged. In this conceptual vacuum, 
models based on simple singlet scalar fields have held sway [T] . The simplest of these are chaotic 
inflation models based on exponential or power-law potentials, of which 4 and (f 2 are the only 
renormalizable examples. The WMAP collaboration has made so bold as to claim that such a 
4 model is excluded at the 3-cr level 1 , a conclusion which would merit further support [31 E]- 
Nevertheless, it is clear that a 2 model would be favoured. 

We reconsider in this paper the possibility that the inflaton could in fact be related to the 
other dramatic recent development in fundamental physics, namely the discovery of neutrino 
masses j^]. The simplest models of neutrino masses invoke heavy singlet neutrinos that give 
masses to the light neutrinos via the seesaw mechanism [B]. The heavy singlet neutrinos are 
usually postulated to weigh 10 10 to 10 15 GeV, embracing the range where the inflaton mass 
should lie, according to WMAP et al. In supersymmetric models, the heavy singlet neutrinos 
have scalar partners with similar masses, sneutrinos, whose properties are ideal for playing the 
inflaton role 7J. In this paper, we discuss the simplest scenario in which the lightest heavy 
singlet sneutrino drives inflation. This scenario constrains in interesting ways many of the 18 
parameters of the minimal seesaw model for generating three non-zero light neutrino masses. 

This minimal sneutrino inflationary scenario (i) yields a simple \m 2 4> 2 potential with no quartic 
terms, with (ii) masses m lying naturally in the inflationary ballpark. The resulting (iii) spectral 
index n s , (iv) the running of n s and (v) the relative tensor strength r are all compatible with the 
data from WMAP and other experiments [2]. Moreover, fixing m ~ 2 x 10 13 GeV as required by 
the observed density perturbations (vi) is compatible with a low reheating temperature of the 
Universe that evades the gravitino problem [H], (vii) realizes leptogenesis [Hlllllj in a calculable 
and viable way, (viii) constrains neutrino model parameters, and (ix) makes testable predictions 
for the flavour-violating decays of charged leptons. 

The main features of our scenario are the following. First, reheating of the Universe is now due 
to the neutrino Yukawa couplings, and therefore can be related to light neutrino masses and 
mixings. Secondly, the lepton asymmetry is created in direct sneutrino-inflaton decays [TU] . 
There is only one parameter describing the efficiency of leptogenesis in this minimal sneutrino 
inflationary scenario in all leptogenesis regimes - the reheating temperature of the Universe - to 

1 This argument applies a fortiori to models with </>™ >4 potentials. 
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which the other relevant parameters can be related. This should be compared with the general 
thermal leptogenesis case [3 El El E3] which has two additional independent parameters, 
namely the lightest heavy neutrino mass and width. Thirdly, imposing the requirement of 
successful leptogenesis, we calculate branching ratios for \i — > and r — > /ry [Hj, and the 
CP- violating observables ^H] like the electric dipole moments of the electron and muon 16J. All 
these leptonic observables, as well as leptogenesis, are related to the measured neutrino masses 
via a parametrization with a random orthogonal matrix ^7j. We show that, in the minimal 
scenario discussed here, successful leptogenesis implies a prediction for \i — > cy in a very narrow 
band within about one order of magnitude of the present experimental bound, whilst r — > 7x7 
might be somewhat further away. 

Other sneutrino inflationary scenarios could be considered. For example, the inflaton might be 
one of the heavier singlet sneutrinos, or two or more sneutrinos might contribute to inflation, 
or one might play a role as a curvaton ^H]- These alternatives certainly merit consideration, 
though they would in general be less predictive. We find it remarkable that the simplest 
sneutrino inflationary scenario considered here works as well as it does. 



2 Chaotic Sneutrino Inflation 

We start by reviewing chaotic inflation pQ with aV = |m 2 2 potential - the form expected for 
a heavy singlet sneutrino - in light of WMAP 0. Defining M P = l/^8nG N ~ 2.4 x 10 18 GeV, 
the conventional slow-roll inflationary parameters are 

l.^/^'V 2M 2 w2 (V"\ 2M 2 P , , (VV"'\ 

where (pi denotes the a priori unknown inflaton field value during inflation at a typical CMB 
scale k. The overall scale of the inflationary potential is normalized by the WMAP data on 
density fluctuations: 

A2r= 2^M 2 1 = 2 - 95x1 °~ 9A : A = - 77±0 - 07 ' ( 2 ) 



y\ = MpVe x 24vr 2 x 2.27 x 1Q- 9 = 0.027M P x ei, (3) 



yielding 

corresponding to 

m f <))! = 0.038 x M| (4) 

in any simple chaotic <p 2 inflationary model, such as the sneutrino model explore here. The 
number of e-foldings after the generation of the CMB density fluctuations observed by COBE 
is estimated to be 
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where prh is the energy density of the Universe when it is reheated after inflation. The second 
term in (jSj) is negligible in our model, whereas the third term could be as large as (—8) for 
a reheating temperature Trh as low as 10 6 GeV. Conservatively, we take iV ~ 50. In a 2 
inflationary model, this implies 

1 ^ 2 



N = 1172* 50 > ( 6 ) 



4MJ, 



corresponding to 



<j>] ~ 200 x M 2 P . (7) 

Inserting this requirement into the WMAP normalization condition ©, we find the following 
required mass for any quadratic inflaton: 

m ~ 1.8 X 10 13 GeV. (8) 

As already mentioned, this is comfortably within the range of heavy singlet (s)neutrino masses 
usually considered, namely m N ~ 10 10 to 10 15 GeV. 

Is this simple <fi 2 sneutrino model compatible with the WMAP data? The primary CMB 
observables are the spectral index 

8M 2 

n s = l-6e + 2r] = l ~ 0.96, 9 

m 



the tensor-to scalar ratio 



r = ^ = 166 = ^^0.16, (10) 



and the spectral-index running 



dn x 2 r, . oi _ . 32M 4 



d\nk 3 



n s - \Y - V + 2£ = ~8x 10" 4 . (11) 



The value of n s extracted from WMAP data depends whether, for example, one combines them 
with other CMB and/or large-scale structure data. However, the (ft 2 sneutrino model value 
n s ~ 0.96 appears to be compatible with the data at the \-o level. The <fi 2 sneutrino model 
value r ~ 0.16 for the relative tensor strength is also compatible with the WMAP data. One 
of the most interesting features of the WMAP analysis is the possibility that dn s /dlnk might 
differ from zero. The <fi 2 sneutrino model value dn s /d\nk ~ 8 x 10~ 4 derived above is negligible 
compared with the WMAP preferred value and its uncertainties. However, dn s /d\nk = 
appears to be compatible with the WMAP analysis at the 2-a level or better, so we do not 
regard this as a death-knell for the (f 2 sneutrino model 2 . 



2 In fact, we note that the favoured individual values for n s ,r and dn s /d\nk reported in an independent 
analysis 0] all coincide with the </> 2 sneutrino model values, within the latter's errors! 
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3 Reheating and Leptogenesis 



Before addressing leptogenesis in this sneutrino model for inflation in all calculational details, 
we first comment on the reheating temperature T R h following the inflationary epoch. Assuming, 
as usual, that the sneutrino inflaton decays when the the Hubble expansion rate H ~ m, and 
that the expansion rate of the Universe is then dominated effectively by non-relativistic matter 
until H ~ T^, where is the inflaton decay width, we estimate 



/ 90 \ 4 / 

7W=(^) y/W, (12) 

where is the number of effective relativistic degrees of freedom in the reheated Universe. In 
the minimal sneutrino inflation scenario considered here we have <fi = N\, m = and 

T^T Nl = ^-(Y u YJ) 11 M Nl , (13) 

where Y u is the neutrino Dirac Yukawa matrix. If the relevant neutrino Yukawa coupling 
(Y v Y}) n ~ 1, the previous choice m = M Nl ~ 2 x 10 13 GeV would yield T RH > 10 14 GeV, 
considerably greater than m itself 3 . Such a large value of T R h would be very problematic for 
the thermal production of gravitinos [8 . However, it is certainly possible that (Y u Y^)u <C 1, 
in which case T R h could be much lower, as we discuss in more detail below. Alternatively, 
one may consider more complicated scenarios, in which three sneutrino species may share the 
inflaton and/or curvaton roles between them. 

We now present more details of reheating and leptogenesis. In general, inflaton decay and the 
reheating of the Universe are described by the following set of Boltzmann equations ^Hj 



dt 
dp R 
dt 



-AHp R + I> , (14) 

jp 

H = m = ^/8nG N (p^ + p R )/3, (15) 



where p^ is the energy density of the inflaton field, p R describes the energy density of the 
thermalized decay products and essentially defines the temperature via 

7T 2 

PR = ^9*T 4 , (16) 

H is the Hubble constant and GV is the Newton constant. Thus reheating can be described by 
two parameters, the reheating temperature ()12J) . which is the highest temperature of thermal 

3 Even such a large value of (Y^Y^n would not alter significantly the <f> 2 sneutrino model prediction for 
dn s / d\nk. 
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plasma immediately after reheating is completed, and the initial energy density of the inflaton 
field 



2„ rpS 



irg 



(17) 



which determines the maximal plasma temperature in the beginning of the reheating process. 
In the following we use the parameter 



M, 



T 



to parametrize temperature. 

The set of Boltzmann equations describing the inflaton decay and reheating, the creation and 
decays of thermal heavy neutrinos and sneutrinos, and the generation of a lepton asymmetry, 
is given by 



HZz 



Z d7 
dY Nl 



3p4> ~^<t>P<j) 



dz 



3F<j>p<j> 
'^P~ 



R 



HZz- 



dY, 



JV 4 



dz 



HZz- 



dY, 



N- 



dz 



zH ' 

Y/Vi (remaining) 



-Yft (remaining) 

4:pR + S 

- - Yfi_ (remaining) 



HZz 
HZz 



dY^ 

dz 
dY Ls 
dz 

H 



sr<t>P<t> 

' ^PR 

'^P~ 



Y Lf + 



Y Ls + 



2sM N ; €l ~ l( remamm 9) 



R 



■e\ (remaining) 

2sMtv, s 



where 



MG N (p (j} + p R )/3, 



Z = 1 - 



AHp 



R 



(19) 
(20) 

(21) 

(22) 

(23) 

(24) 
(25) 

(26) 



N± = Ni±Nl, and Y Nl , Y^ ± , Y Lf , Y Ls , denote the number-density-to-entropy ratios, Y = n/s, 
for the heavy neutrinos, sneutrinos and lepton asymmetries in fermions and scalars, respectively. 
The terms denoted by remaining are the usual ones for thermal leptogenesis, and can be 
obtained from [TJ] by using H(Mn 1 ) = z 2 H. We do not write out their lengthy expressions in 
full here. The first terms on the r.h.s. of (|2()H24j) are the dilution factors of Y = n/s due to 
entropy production in the inflaton = N x decays described by (fT§j). The second terms on the 
r.h.s. of (|23|). ()24|) describe lepton asymmetry generation in the decays of the coherent inflaton 
field. Identifying and studying the parameter space in which leptogenesis is predominantly 
direct is one of the aims of this paper. 
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Figure 1: Lower bound (solid curve) on as a function offh\ forYs > 7.8 x 10 _11 ; assuming 
a maximal CP asymmetry e™ ax (M Nl ). Successful leptogenesis is possible in the area above the 
solid curve. In the area bounded by the red dashed curve, leptogenesis is entirely thermal. 



We are now ready to study (|19H25jl . First we work out general results on reheating and lepto- 
genesis in the sneutrino inflation scenario, allowing to vary as a free parameter. In this 
case, the reheating and leptogenesis efficiency is described by two parameters, namely M Nl and 
a parameter describing the decays of the sneutrino inflaton. This can be chosen to be either 
rhi = (XvY^)iiv 2 sin 2 or, more appropriately for this scenario, the reheating temperature 

of the Universe Trh given by (|T2"|) . For the CP asymmetry in (s)neutrino decays, we take the 
maximal value for hierarchical light neutrinos, given by (2Hj: 



This choice allows us to study the minimal values for Mn x and Trh allowed by leptogenesis. 
Later, we will focus our attention on exact values of e\ |2~T] . 

Solutions to (J19H25|) are presented in Figs. [T] and El We plot in Fig. [T] the parameter space 
in the (M^, mi) plane that leads to successful leptogenesis. This parameter space has three 
distinctive parts with very different physics. 

In the area bounded by the red dashed curve, denoted by A, leptogenesis is entirely thermal. 
This region has been studied in detail in Whatever lepton asymmetry is generated initially 
in the decay of the sneutrino inflaton is washed out by thermal effects, and the observed baryon 
asymmetry is generated by the out-of-equilibrium decays of thermally created singlet neutrinos 





(27) 
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Figure 2: The solid curve bounds the region allowed for leptogenesis in the (Trh, M^) plane, 
again obtained assuming Yb > 7.8 x 10 -11 and the maximal CP asymmetry e r p ax (M^ 1 ). In the 
area bounded by the red dashed curve leptogenesis is entirely thermal. 



and sneutrinos. As seen in Fig. |21 in our scenario this parameter space corresponds to high 
M Nl and high T RH values. 

The area B below the dashed curve and extending down to the minimum value = 4 x 10 6 
GeV in Fig. ^ is the region of parameter space where there is a delicate cancellation between 
direct lepton asymmetry production in sneutrino inflaton decays and thermal washout. This 
region cannot be studied without solving the Boltzmann equations numerically. However, it 
roughly corresponds to Trh ~ as seen in Fig. El 

The area denoted by C has Trh <C M^ 1 . Since the maximal CP asymmetry scales with M^, 
the line presented corresponds to a constant reheating temperature. Notice that in Fig. ^ this 
line is terminated at fh\ = 10~ 7 . As seen in Fig. it continues linearly to high values of 
Mn 1 . In this area, leptogenesis is entirely given by the decays of cold sneutrino inflatons, a 
scenario studied previously in jlOj . In this case the details of reheating are not important for 
our analyses. To calculate the lepton asymmetry to entropy density ratio Yl — Ul/s in inflaton 
decays we need to know the produced entropy density 



and to take into account that inflaton dominates the Universe. In this case one obtains ^0] 



s = 




45 



(28) 




(29) 
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where €\ is the CP asymmetry in <p = N\ decays. The observed baryon asymmetry of the 
Universe gives a lower bound on the reheating temperature Trh > 10 6 GeV. 

We consider now the most constrained scenario in which the inflaton is the lightest sneutrino, 
which requires Mn 3 > M^ 2 > — 2 x 10 13 GeV. This implies that our problem is completely 
characterized by only one parameter, either rh\ or Trh- As we see in both Figs. ^and|21 the 
line for ~ 2 x 10 13 GeV traverses both the regions A, and C, the former corresponding to 
high Trh, as seen in Fig. 12 However, Trh may also be low even in the minimal seesaw model, 
as seen in Fig. El 

The cosmological gravitino problem suggests that Trh ^ 10 8 GeV might be the most interest- 
ing, which would correspond to very small rh\, far away from the thermal region A and deep in 
the region C where leptogenesis arises from the direct decays of cold sneutrinos. We concentrate 
on this option here. This limit requires very small Yukawa couplings (Y v YJ)n ^ 10 -12 , whilst 
other Yukawa couplings can be 0(1). This possibility may be made natural, e.g., by postulating 
a Z 2 matter parity under which only N% is odd. In this case, the relevant Yukawa couplings (Yv)) 
all vanish, but a Majorana mass for N\ is still allowed. A more sophisticated model postulates 
a Z7 discrete family symmetry with charges Yrn = (4, 0, 0) for the Ni, (2, 1, 1) for the 5 repre- 
sentations of SU(5), and (2, 1, 0) for the 10 representations of SU(5). Assuming a gauge-singlet 
field $ with Y FN = -1 and ($) = e, we find M { = 0(e, 1, 1) and %)) = C(e 6 , e 5 , e 5 ), whilst the 
other Yukawa couplings are 0{1), 0(e) or 0(e 2 ). If e ~ 1/17, the (Y v )) are sufficiently small 
for our purposes, whilst the quark and lepton mass matrices are of desirable form. Doubtless, 
one could construct better models with more effort, but this example serves as an existence 
proof for a low value of Trh in our scenario. 



4 Leptogenesis Predictions for Lepton Flavour Violation 

In this Section, we relate the results of the previous section on direct leptogenesis to light 
neutrino masses, and make predictions on the lepton-flavour-violating (LFV) decays. Thermal 
leptogenesis in this context has been extensively studied recently [22123 123] • We first calculate 
neutrino Yukawa couplings using the parametrization in terms of the light and heavy neutrino 
masses, mixings and the orthogonal parameter matrix given in ^7j. This allows us to calculate 
exactly the baryon asymmetry of the Universe, since we know the CP asymmetry e\ and 
the reheating temperature of the Universe Trh- For neutrino parameters yielding successful 
leptogenesis, we calculate the branching ratios of LFV decays. 

There are 18 free parameters in the minimal seesaw model with three non-zero light neutrinos, 
which we treat as follows. In making Fig. El we have taken the values of #12, #23, Am^ and 
Am 2 ^ from neutrino oscillation experiments. We randomly generate the lightest neutrino mass 
in the range < m x < 0.01 eV and values of 9 13 in the range < #13 < 0.1 allowed by 
the Chooz experiment [2H], as we discuss later in more detail. Motivated by our previous 
discussion of chaotic sneutrino inflation, we fix the lightest heavy singlet sneutrino mass to be 
Mi = 2 x 10 13 GeV, and choose the following values of the heavier singlet sneutrino masses: 
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M 2 = 10 14 GeV or M 2 = 5 x 10 14 GeV, and M 3 in the range 5 x 10 14 to 5 x 10 15 GeV, as we 
also discuss later in more detail. This accounts for nine of the 18 seesaw parameters. 

The remaining 9 parameters are all generated randomly. These include the three light-neutrino 
phases - the Maki-Nakagawa-Sakata oscillation phase and the two Majorana phases. Specifi- 
cation of the neutrino Yukawa coupling matrix requires three more mixing angles and three 
more CP- violating phases that are relevant to leptogenesis, in principle. The plots in Fig. El are 
made by sampling randomly these nine parameters. We apply one constraint, namely that the 
generated baryon density falls within the 3 — o range required by cosmological measurements, 
of which the most precise is now that by WMAP: 7.8 x KT 11 < Y B < 1.0 x 1CT 10 0. 

Making predictions for LFV decays also requires some hypotheses on the parameters of the 
MSSM. We assume that the soft supersymmetry-breaking mass parameters mo of the squarks 
and sleptons are universal, and likewise the gaugino masses m^, and we set the trilinear soft 
supersymmetry-breaking parameter A Q = at the GUT scale. Motivated by g M — 2, we assume 
that the higgsino mixing parameter \x > 0, and choose the representative value tan/? = 10. 
We take into account laboratory and cosmological constraints on the MSSM, including limits 
on the relic density of cold dark matter. WMAP provides the most stringent bound on the 
latter, which we assume to be dominated by the lightest neutralino \: 0.094 < VL x h 2 < 0.129. 
For tan/5 = 10, the allowed domain of the (mi/2,mo) plane is an almost linear strip extending 
from (mi/2,m ) = (300,70) GeV to (900,200) GeV [26 . For illustrative purposes, we choose 
(mift, mo) = (800, 170) GeV and comment later on the variation with mu2- 

Panel (a) of Fig. El presents results on the branching ratio BR for fi — > ej decay. We see 
immediately that values of Trh anywhere between 2 x 10 6 GeV and 10 12 GeV are attainable in 
principle. The lower bound is due to the lower bound on the CP asymmetry, while the upper 
bound comes from the gravitino problem. The black points in panel (a) correspond to the 
choice sin #13 = 0.0, M 2 = 10 14 GeV, and 5 x 10 14 GeV < M 3 < 5 x 10 15 GeV. The red points 
correspond to sin# 13 = 0.0, M 2 = 5 x 10 14 GeV, and M 3 = 5 x 10 15 GeV, while the green points 
correspond to sin6 l 13 = 0.1, M 2 = 10 14 GeV, and M 3 = 5 x 10 14 GeV. We see a very striking 
narrow, densely populated bands for BR(/i — > ej), with some outlying points at both larger 
and smaller values of BR(/i — > ej). The width of the black band is due to variation of Mjv 3 
showing that BR(/x — > ej) is not very sensitive to it. However, BR(/i — > ej) strongly depends 
on M N2 and sin # 13 as seen by the red and green points, respectively. Since BR(/i — > ej) scales 
approximately as m^^ 2 , the lower strip for sin# 13 = would move up close to the experimental 
limit if mi/2 ~ 500 GeV, and the upper strip for sin#i 3 = 0.1 would be excluded by experiment. 

Panel (b) of Fig. El presents the corresponding results for BR(r — * /ry) with the same colour 
code for the parameters. This figure shows that BR(r — > /ry) depends strongly on M N3 , while 
the dependence on sin# 13 and on M^ 2 is negligible. The numerical values of BR(r — > /ry) 
are somewhat below the present experimental upper limit BR(r — > /ry) ~ 10~ 7 , but we note 
that the results would all be increased by an order of magnitude if mu^ ~ 500 GeV. In this 
case, panel (a) of Fig. El tells us that the experimental bound on BR(yU — > ej) would enforce 
sin^is 0.1, but this would still be compatible with BR(r — > /ry) > 10~ 8 . 

As a result, Fig. El strongly suggests that fixing the observed baryon asymmetry of the Universe 
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Figure 3: Calculations of BR(fi — > e7) and BR(r — » /ry) on Ze/£ and ng/rf panels, respectively. 
Black points correspond to sin #13 = 0.0, M 2 = 10 14 GeV, and 5 x 10 14 CeF < M 3 < 5 x 
10 15 GeK iZed points correspond to sinfl 13 = 0.0, M 2 = 5 x 10 14 GeV, and M 3 = 5 x 10 15 GeF, 
u>/w/e green points correspond to sin# 13 = 0.1, M 2 = 10 14 GeV, and M 3 = 5 x 10 14 GeV. 



for the direct sneutrino leptogenesis (Tbh < 2 x 10 12 GeV < MjvJ implies a prediction for the 
LFV decays provided M N2 and/or M Na are also fixed. This observation can be understood in 
the case of hierarchical light and heavy neutrino masses. Consider first fi — ► cy for sin^ = 0. 
It turns out that the iV 2 couplings dominate in (Y u Yj) 2 i which determines BR(,u — > ej). Also, 

the M N2 term dominates in t\ which implies Y B ~ (Y u Yj) 2 i/ \f (Y u Yj) u , because cancellations 
among the phases are unnatural. In the parametrization with the orthogonal matrix R, this 
implies Yb ~ -^23/^22- If fine tunings are not allowed, the requirement Trh < fixes 
R23/R22 and therefore relates Yb to fi — > ej. For more general cases, the behaviour of BR(/i — > 
ej) is more complicated and additional contributions occur. However, those new contributions 
tend to enhance BR(/i — > e^) , as exemplified in Fig. |3] by the green dots. 

The behaviour of BR(r — > fi'-f) is simpler. To leading order in the largest parameters, r — > fi'j 
depends on the N 3 couplings and mass, leading to (Y u Yj) 3 2 ~ (^)33^33^23) independently of 
leptogenesis results. 

We have to stress here that such definite predictions for LFV processes can always be avoided 
by fine tuning the neutrino parameters, as seen by several scattered points in Fig. El Points 
with small BR(/i — > e^) can be systematically generated using the parametrization of Y v by a 
Hermitian matrix [21] , and the predictions for the LFV decays thereby washed away. However, 
in this case, the M^. are outputs of the parametrization, and cannot be fixed as required by the 
present analyses of sneutrino inflation. Therefore the parametrization [2Z| is not appropriate 
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for our leptogenesis scenario. Finally, we comment that such fine tunings are impossible in 
simple models of neutrino masses [2%] . 

Another possibility for avoiding the LFV predictions is to allow the heavy neutrinos to be 
partially degenerate in mass, which enhances the CP asymmetries j2H]- in supersymmetric 
models, this possibility was considered in [2Ti] . 

In addition to the quantities shown in Fig. El we have also examined BR(r — > ej), which is 
always far below the present experimental bound BR(r — > erf) ~ 10~ 7 , and the electron and 
muon electric dipole moments. We find that d e < 10~ 33 e cm, in general, putting it beyond the 
foreseeable experimental reach, and \d^/d e \ ~ m M /m e , rendering also unobservably small. 

5 Alternative Scenarios and Conclusions 

We have considered in this paper the simplest sneutrino inflation scenario, in which the inflaton 
4> is identified with the lightest sneutrino, and its decays are directly responsible for leptogenesis. 
We find it remarkable that this simple scenario is not already ruled out, and have noted the 
strong constraints it must satisfy enable it to make strong predictions, both for CMB observables 
and LFV decays. These might soon be found or invalidated. In the latter case the motivation 
to study more complicated sneutrino inflation scenarios would be increased. 

• One possibility is that inflation might have been driven by a different sneutrino, not the 
lightest one. In this case, the lightest sneutrino could in principle be considerably lighter than 
the 2 x 10 13 GeV required for the inflaton. This would seem to make more plausible a low 
reheating temperature, as suggested by the gravitino problem. However, this problem is not 
necessarily a critical issue, as it can already be avoided in the simplest sneutrino inflation 
scenario, as we have seen. On the other hand, if the lightest sneutrino is not the inflaton, 
leptogenesis decouples from inflationary reheating, and predictivity is diminished. 

• A second possibility is that two or more sneutrinos contribute to inflation. In this case, 
the model predictions for the CMB observables and the sneutrino mass would in general be 
changed. 

• A related third possibility is that one or more sneutrinos might function as a curvaton, which 
would also weaken the CMB and sneutrino mass predictions. 

For the moment, we do not see the need to adopt any of these more complicated scenarios, but 
they certainly merit investigation, even ahead of the probable demise of the simplest sneutrino 
inflation scenario investigated here. 
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